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We report a measurement of the 5° lifetime in the semileptonic decay channel 5° -^ D~ fi'^vX 
(and its charge conjugate), using approximately 0.4 fb^^ of data collected with the DO detector 
during 2002-2004. Using 5176 reconstructed D~ n'^ signal events, we have measured the 5° lifetime 
to be t{Bs) = 1.398 ± 0.044 (stat) J^qqjs (syst) ps. This is the most precise measurement of the Bg 



lifetime to date. 
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Measurements of the lifetimes of different b hadrons al- 
low tests of the mechanism of heavy hadron decay. The 
spectator model predicts that all hadrons with the same 
heavy flavor content have identical lifetimes. However, 
observed charm and bottom hadron lifetimes suggest that 
non-spectator effects, such as interference between con- 
tributing amplitudes, are not negligible in heavy hadron 
decays. This implies that a mechanism beyond the sim- 
ple spectator model is required. An effective theory called 
the Heavy Quark Expansion (HQE) 'l'| includes such ef- 
fects and predicts lifetime differences among the different 
bottom hadrons. In particular, a difference of the order 
of 1% is predicted between B^ and J5° mesons. The mea- 
surement of the flavor-specific B'^ lifetime using semilep- 
tonic decays is also useful in determining the decay width 
difference between the light and heavy mass eigenstates 
of the i?" meson, which is an equal mixture of CP eigen- 
states that correspond to mass eigenstates in the absence 
of CP violation in the B'^ system. 

In this Letter, we present a high-statistics measure- 
ment of the -Bg lifetime, using a large sample of semilep- 
tonic B° decays collected in pp collisions at ^/s = 
1.96 TeV with the DO detector at the Fermilab Teva- 
tron Colhder in 2002 - 2004. The data correspond to 
approximately 0.4 fb~^ of integrated luminosity. 5° 
mesons were identified through their semileptonic decay 
S° -^ D^ fi^vX [13|, where the Dj meson decays via 
D^ -^ (j>Tr^ , followed by (/) — > K~^K~. 

The DO detector is described in detail elsewhere [3. 
The detector components most important to this analy- 
sis are the central tracking and muon systems. The DO 
central-tracking system consists of a silicon microstrip 
tracker (SMT) and a central fiber tracker (CFT), both 
located within a 2 T superconducting solenoidal mag- 
net, with designs optimized for tracking and vertexing 
at pseudorapidities \?]\ < 3 and \t]\ < 2.5, respectively 
(where 77 = — ln[tan(0/2)]). A liquid-argon and uranium 
calorimeter has a central section covering pseudorapidi- 
ties up to « 1.1, and two end calorimeters that extend 
the coverage to I77I « 4.2 Q. The muon system is located 
outside the calorimeters and has pseudorapidity cover- 
age I77I < 2. It consists of a layer of tracking detectors 
and scintillation trigger counters in front of 1.8 T toroids, 
followed by two similar layers after the toroids [J|. 

Events with semileptonic B-meson decays were se- 
lected using inclusive single-muon triggers in a three-level 
trigger system. The triggers used did not impose any im- 
pact parameter criterion and were shown to not bias the 
lifetime measurement. Off-line, muons were identified by 
extrapolation of the muon track segments, formed by the 
hits in the muon system, to the tracks found in the cen- 
tral tracking system. Each muon was required to have 



a momentum p > 3 GeV/c and a transverse momentum 
PT > 2 GeV/c. 

The primary vertex of each pp interaction was defined 
by all available well-reconstructed tracks 5!| and con- 
strained by the mean beam-spot position. The latter 
was updated every few hours. The resolution of the re- 
constructed primary vertex was typically 20 /im in the 
transverse plane and 40 /xm in the beam direction. 

To reconstruct D^ -^ (f)Tr^ decays, tracks with pr > 
1.0 GeV/c were assigned the kaon mass and oppositely 
charged pairs were combined to form a candidate. Each 
(j) candidate was required to have a mass in the range 
1.008 - 1.032 GeV/c^, compatible with the reconstructed 
(j) mass at DO. The (j) candidate was then combined with 
another track of pr > 0.7 GeV/c. For the "right-sign" 
combinations, we required the charge of the track to be 
opposite to that of the muon and assigned the pion mass 
to this track. All selected tracks were required to have 
at least one SMT hit and one CFT hit. The three tracks 
selected were combined to form a common vertex (the 
D~ vertex) with a confidence level greater than 0.1%. 
The Dj candidate was required to have pr > 3.5 GeV/c. 

The secondary vertex, where the i?^ decays to a muon 
and a Dj meson, was obtained by finding the intersec- 
tion of the trajectory of the muon track and the fiight 
path of the D~ candidate. The confidence level of that 
vertex had to be greater that 0.01%. The reconstructed 
Dj decay vertex was required to be displaced from the 
primary vertex in the direction of the Dj^ momentum. 

The helicity angle, $, defined as the angle between the 
directions of the K~ and D~ in the (j) rest frame, has a 
distribution proportional to cos^ $. A cut of |cos$| > 
0.4 was applied to further reduce combinatorial back- 
ground, which was found to have a flat distribution. 
In order to suppress the physics background originating 
from £)(*)£)(*) processes [IJ], we required that the trans- 
verse momentum of the muon with respect to the Dj me- 
son, pTrei, exceed 2 GeV/c. The Djfj,~^ invariant mass 
was also restricted to 3.4 — 5.0 GeV/c^, to be consistent 
with a _B-meson candidate. Since the number of tracks 
near the -B^ candidate tends to be small, we required the 
isolation X = p*°*(/Lt+D7)/(p*°*(/Lt+ 1)7 ) + X;pf*) > 0.65, 
where the sum X]Pi°* ^^^ taken over all charged parti- 
cles in the cone •\/(A0)2 + (A*?)^ < 0.5, with Acj) and Arj 
being the azimuthal angle and the pseudorapidity with 
respect to the {iX^D~) direction. The muon, kaon, and 
pion tracks were not included in the sum. 

The lifetime of the S^, r, is related to the decay 
length in the transverse plane, L^y, by L^y = crpT/m, 
where pT is the transverse momentum of the B'^ and 
m is its invariant mass. L^y is defined as the displace- 
ment of the i3g vertex from the primary vertex pro- 
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FIG. 1: The mass distribution of <j!>7r~ candidates. Points 
with errors bars show the "right-sign" D~ fi'^ combinations, 
and the open squares show the corresponding "wrong-sign" 
distribution. The dashed curve represents the result of the fit 
to the "right-sign" combinations. The two peaks are associ- 
ated with the D~ and DJ mesons, respectively. 



jected onto the transverse momentuni of the Djfi^ sys- 
tem. Since the B° meson is not fully reconstructed, 
Pt{B'^) is estimated by pt{D~ ii'^)/K, where the cor- 
rection factor, K = pt{D~ fji^) / pt{B^) is determined 
using Monte Carlo (MC) methods. The quantity used to 
extract the B^ lifetime is called the pseudo-proper decay 
length (PPDL). The correction factor K was applied sta- 
tistically when extracting ct{B^) from the PPDL in the 
lifetime fit. 

In the cases with more than one 5° candidate per 
event, we chose the one with the highest vertex confidence 
level. We also required the PPDL uncertainty to be less 
than 500 /im. The resulting invariant mass distribution 
of the Z?7 candidates is shown in Fig. [TJ The distribu- 
tion for "right-sign" Dj ji^ candidates was fitted using a 
Gaussian, to describe the signal, and a second-order poly- 
nomial, to describe the combinatorial background. A sec- 
ond Gaussian was included for the Cabibbo-suppressed 
D~ -^ 07r~ decay. The best fit result is shown in the 
same figure. The fit yields a signal of 5176 ± 242 (stat) 
± 314 (syst) D~ candidates and a mass of 1958.8 ± 0.9 
MeV/c2. The width of the Dj Gaussian is 22.6 ± 1.0 
MeV/c^. The systematic uncertainty comes from the fit. 
For the D~ meson, the fit yields 1551 events. Figure [1] 
also shows the invariant mass distribution of the "wrong- 
sign" candidates. The observed shift in the D7 mass is 
consistent with known issues associated with the calibra- 
tion of the DO track momenta. The contribution to the 
mass region from reflected states was found to be negli- 
gible. Studies confirmed that this mass shift introduces 
no significant residual bias in the lifetime determination. 

MC samples were generated using Pythia [6] for the 
production and hadronization phase, and EvtGen [7| 
for decaying the b and c hadrons. Branching ratios from 



the PDG have been used when available. Detector accep- 
tance and smearing were taken into account using the full 
DO detector simulation based on GEANT l8|. Generated 
MC signal samples include contributions from D^fi^v, 
D*-iJL+v, D^-fi+u, D'-fi+v, and Di*^^T+iy. 

Apart from the background due to combinatorial pro- 
cesses such as a prompt muon and an identified Dj me- 
son, there are several real physics processes that pro- 
duce a muon and a Dj meson, where neither comes 
from the semileptonic decay of the _B^ meson. These 
"right-sign" D~ fi^ combinations are included in the sig- 
nal sample and are defined as "physics backgrounds." 
Prompt D'^ mesons from cc production at the interac- 
tion point can combine with high-p^ muons generated 
either via direct production or in charm decays. These 
cc background events are expected to have very short 
lifetimes and thus could introduce a significant bias in 
the i?° lifetime measurement. Backgrounds that origi- 
nate from B mesons and provide the D^ 11^ final state, 
but not via the semileptonic decay B^ -^ D^ fi^vX, are 
called non-_B^ backgrounds. This kind of background is 
expected to have a relatively long lifetime, thus its ef- 
fect on the _B" lifetime fit is smaller than that of the 
charm background. There are three sources of such 
events: B" -> oi*'^^ D^*^+X, B- -^ oi*'^^ D^*'>"X, and 
B^ -^ DsD'^*^X, where the charm meson accompa- 
nying the Dg , which decays to (/>7r~, decays semilep- 
tonically. The momentum of the muon coming from the 
decay of the D'-*^ is softer than that for the signal, since 
it comes from the decay of a secondary charm hadron. 
This implies that the contribution of these modes to the 
signal sample is reduced by the kinematic cuts. We found 
the fractional contribution of the backgrounds to the sig- 
nal region to be (10.0 ± 7.0)% for cc background and 
(11.3^3:6)% for non-B^ backgrounds. 

The lifetime of the B^ was found using a fit to the 
PPDL distribution. We defined a signal sample using the 
-D7 mass distribution in the region from 1913.6 MeV/c^ 
to 2004.0 MeV/c^, corresponding to ±2ct from the fit- 
ted mean mass. The PPDL distribution of the combi- 
natorial background events contained in the signal sam- 
ple was defined using "right-sign" events from the D~ 
sidebands (1755.3 - 1800.5 MeV/c^ and 2117.1 - 2162.3 
yieV /(?) and "wrong-sign" events between 1755.3 and 
2162.3 MeV/c^. The combinatorial background due to 
random track combinations was modeled by the sideband 
sample events. This assumption is supported by the mass 
distribution of the "wrong-sign" combinations where no 
enhancement is visible in the D~ mass region. 

The PPDL distribution obtained from the signal sam- 
ple was fitted using an unbinned maximum log-likelihood 
method. Both the B'^ lifetime and the background shape 
were determined in a simultaneous fit to the signal and 
background samples. The likelihood function C is given 
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where Ns, Nb are the number of events in the signal 
and background samples and fsig is the ratio of Dj sig- 
nal events obtained from the Dj mass distribution fit to 
the total number of events in the signal sample. To con- 
strain fsig, we factored in an additional likelihood term 
using the number of Dj signal events observed from the 
invariant mass distribution, and its uncertainty, Csig- 

Since the current world average width difference be- 
tween the light and heavy mass eigenstates (AF^) of the 
B'^ system is small [10[ compared with the current preci- 
sion of the data, we used for the signal probability distri- 
bution function (PDF), J-]j„, a normalized single expo- 
nential decay function convoluted with a Gaussian reso- 
lution function. The if-factor correction was also convo- 
luted with the exponential decay function. Since a priori, 
we do not know the decay length uncertainty, which we 
estimated on an event-by-event basis, an overall global 
scale factor, s, was introduced as a free parameter in the 
B'^ lifetime fit. The events from non-i?° background were 
taken into account in the fit by including similar PDFs 
to those in the signal but using fixed parameters accord- 
ing to the world-average values [9| . A different iiT-factor 
distribution was also used for each process. To take into 
account the cc background, we used a Gaussian distri- 
bution with fixed parameters. These contributions were 
evaluated and parametrized using MC methods following 
a similar procedure as for the signal evaluation. 

The combinatorial background sample, ^l^-j^, was para- 
metrized using a Gaussian distribution function for the 
resolution plus several exponential decays: two for the 
negative values in the PPDL distribution (one short and 
one long component) and two for the positive values of 
the distribution. 

Figure [2] shows the PPDL distribution of the D^ fi^ 
signal sample with the fit result superimposed (dashed 
curve). The dotted curve represents the sum of the back- 
ground probability function over the events in the signal 
sample. The B'^ signal is represented by the filled area. 

To test the resolutions, pulls, fitting, and selection cri- 
teria, we performed detailed studies using MC samples 
and found no significant bias in our analysis procedure. 
In order to study the stability of the B'^ lifetime measure- 
ment, we split the data sample into two parts according 
to different kinematic and geometric parameters, com- 
pared the fitted results, and found the lifetimes consis- 
tent within their uncertainties. We also varied the selec- 
tion criteria and mass fit ranges, and did not observe any 
significant shifts. We performed an extensive study of 
our fitting procedure, looking for any possible bias using 
MC ensembles with statistics of the size of our dataset 
and distributions as those in data. These samples were 
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FIG. 2; Pseudo-proper decay length distribution for D~ jj,'^ 
candidates with the result of the fit superimposed as the 
dashed curve. The dotted curve represents the combinato- 
rial background and the filled area represents the B^ signal. 



fitted, and the mean and width of the distributions of ex- 
tracted parameters were found to be consistent with the 
fits to data. One final check of the procedure involved 
performing a similar lifetime fit to a control sample de- 
fined by the Cabibbo-suppressed decay D^ -^ ^tt^, (see 
Fig. [T|). We found that 89.1% of the sample comes from 
B° -^ D-^i+X, and the B° lifetime to be 1.541 ± 0.093 
ps, where the uncertainty is statistical only. This re- 
sult is ill good agreement with the world average i3° life- 
time [i,[3- 

We considered and evaluated various sources of system- 
atic uncertainties. The major contributions come from 
the determination of the combinatorial background, the 
model for the resolution, and the physics background. To 
determine the systematics due to the uncertainty on the 
combinatorial background, we tested other assumptions 
on the background samples: we used just the events in 
the sidebands, just the events in the wrong-sign combi- 
nations, and removed either the right sideband or the 
left sideband samples. We also modified the definitions 
of those samples, changing the mass window sizes and 
positions. The largest difference in ct observed in these 
variations of background modeling was 4.3 /um, which 
was taken to be the systematic uncertainty due to this 
source. The effect of uncertainty in the resolution of the 
decay length was studied using an alternative global scale 
factor, s. We repeated the lifetime fit with fixed values 
of s obtained from MC samples and from a different life- 
time analysis 



ll|. Using a variation of the resolution 



scale by a factor of two beyond these bounds, we found a 
3.7 /im variation in ct. The uncertainty from the physics 
background was evaluated by varying the branching frac- 
tions of the different processes as well as the shapes of 
the lifetime templates, as given by their known lifetime 
values [9]. The variations were within one standard de- 
viation in each case. Assuming no correlation between 



them, we added the effects of all the variations in quadra- 
ture and found a total contribution of _4'2 A*™.- Using a 
similar procedure, we evaluated the uncertainty coming 
from the determination of the cc background and found 
a difference of to'g M™- 

To evaluate the uncertainty associated with the K fac- 
tor determination, we modified the kinematics of the 
event using a different decay model, a different pt spec- 
trum for the b quark, and a different pr spectrum for the 
muon. We also varied the amount of each component, 
according to their uncertainty, of the B^ -^ D~^~^X sig- 
nal. In each case, the K factor was re-evaluated and the 
fit repeated. We added all K factor variation effects in 
quadrature and found a total uncertainty of I2 i M™- 

There are two requirements in our selection method 
that could potentially change the final result by altering 
the shape of the PPDL distribution: pTrei > 2 GeV/c 
and the positive displacement from the primary vertex 
of the reconstructed Dj decay vertex. Using MC meth- 
ods, we evaluated their effects by removing them one at 
a time. The largest variation observed was _q'3 fim, and 
the selection efficiency is flat as a function of proper decay 
time. The effect of a possible misalignment of the SMT 
system was tested in Ref. [11]. We repeated the study 
using MC signal samples and observed the same shift of 
CT = 2 /im, which was taken as a systematic uncertainty 
due to possible misalignment. The total systematic un- 
certainty from all of these sources added in quadrature is 

In summary, using an integrated luminosity of ap- 
proximately 0.4 fb^^, we have measured the i?" life- 
time in the decay channel D~^~^i/X to be t{B'^) = 



1.398 ±0.044 (stat) toml (syst) ps. 



Note that this mea- 
surement takes Ar^ equal to zero. The extraction of the 
average lifetime fs for AF^ 7^ is straight forward [lO|. 
The result is in good agreement with previous experi- 
ments as well as the current world average value for all 
flavor-specific decays, t{B°^) = 1.442 ± 0.066 ps 0,111. 
Our B° lifetime measurement is the most precise to date 
and exceeds the precision of the current world average 
measurement t{B^)pdg — 1.461 ± 0.057 ps ^], where 
semileptonic and hadronic decays were combined. This 
measurement is approximately 2.5cr away from the B^ 
lifetime, more than the 1% predicted by HQE. 
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